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■ The three-dimensional stability problem of a stretched stationary vortex is addressed in this 

^-H ' letter. More specifically, we prove that the discrete part of the temporal spectrum is only associated 

with two-dimensional perturbations. 
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I. INTRODUCTION 



Numerical simulations as well as real experiments Qj indicate that vorticity in turbulent flows concentrates in 
' localized regions such as filaments which are fairly well described by stretched vortices such as the celebrated Burgers 
\ vortex solution S or Moffatt, Kida & Ohkitani's asymptotic solution If one agrees that these local structures are 
PsJ 1 important dynamical objects of the global turbulent field, their temporal stability with respect to generic perturbations 
should be addressed. So far, this problem has only been studied for Burgers' vortex and, in that case, for purely 
two-dimensional perturbations. Robinson & Saffman H provided an analytical solution for low-Reynolds numbers 
and their result were later numerically extended by Prochazka & Pullin j^] up to Reynolds numbers Re = 10 4 . 
These papers on Burgers vortex indicate that the temporal spectrum associated with two-dimensional perturbations 
is discrete and corresponds to damped modes. On the contrary, the general stability problem of stretched vortices 
has not been tackled yet. Except for the 2D stability analysis of axisymmetric Burgers vortex, it does not reduce to 
a classical eigenvalue problem with a single ODE to solve. Indeed, infinitesimal 3D perturbations are affected by the 
d ', presence of stretching along the vortex axis which precludes the reduction of the problem by Fourier analysis. 

In this paper, we prove that the discrete part of the temporal spectrum is only associated with two-dimensional 
perturbations. In section |n[ the stability problem is introduced and particular time-dependent solutions are exhibited. 



Their existence imposes conditions on the 3D temporal mode structure. In section III, these conditions are shown to 
be consistent only for modes independent of the vortex axis coordinate. 

II. MODIFIED FOURIER DECOMPOSITION 

Let us consider a stationary velocity field Uo = (Uo, Vq, Wo) of the form 

U = ^(x,y) + U v {x,y) , (2.1a) 

V = ^(x,y)+V v (x,y) , (2.1b) 

W =7Z, (2.1c) 

where (U V ,V V ) and (j^(x,y), ^(x,y),jz) respectively stand for a localized rotational field and a global velocity 
field satisfying V 2 </> = —7. Such an expression represents a stationary stretched vortex aligned with the z-axis and 
subjected to a global strain field. In the sequel, the strain rate 7 along the vortex axis is assumed to be positive. 

The structure of such a solution is governed by the balance betw een s tretching due to the global strain field 
and viscous diffusion. In particular, the core size should scale as y/v/j where v is the kinematic viscosity. In 
the simplest case of an axisymmetric strain (^(x,y), ^(x,y),jz) = (—jx/2, —7^/2,72), one recovers the Burgers 
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solution. Other examples are Robinson & Saffman jjj and Moffatt, Kida & Ohkitani || solutions which correspond to 
t he n on-axisymmetric case at small and large Reynolds numbers respectively. In the subsequent analysis, expression 
( |2.1a -c) is considered as the basic flow where (U Vl V v ) and cf> are not specified. 

The dynamics of pressure and velocity infinitesimal perturbations (u,p) around ( 2.1a -c) is described by the linear 
system : 



d t u + U .Vu + (u.V)Uo = -Vp + f V 2 u . 

V.u = 



(2.2a) 
(2.2b) 



The above system being homogeneous with respect to time, one might look at the temporal spectrum of such 
system. Modes belonging to the discrete part of this spectrum read: 



(u w ,_p w ) = (v ul (x,y,z),q a ,(x,y,z))e 



(2.3) 



Inserting ( |2.3| ) into (2.2a-b), one obtains equations for v w (a;, y, z) and q u (x, y, z) which are non-separable with respect 
to any spatial variable. In such a case, the use of standard Fourier analysis does not simplify any further the problem. 
However, the z-dependence in equation (2.2a) only appears through the uniform strain along the z-axis i.e. through 
the term jzd z . Such a simple dependence allows to search for time-dependent solutions which are different from ( [2.3] ). 
These are "generalized Fourier modes" in the z direction with a time-dependent wavenumber : 



(u,p) = (u(x,y,t,ko),p(x,y,t,k ))e ik ^>e- v f> {s)r ^ , 



(2.4) 



where the initial wavenumber condition kg — fc(0) is a free parameter. Indeed, as soon as the time evolu tion of k(t) 
is appropriately chosen, more precisely if k(t) = fcoe -7 *, the non- homogeneous term in z is removed in ( [2. 2a ). The 
system ( |2.2a ,b) is then reduced to a couple of equations homogeneous in z, which describes the time evolution of the 
modified Fourier modes components u x and u y : 



e 2lt d 

(£ + d x ll )u x + dyU Uy = ^"g~(^ + 2 7) 

e 2 ^ d 

(C + dyV )u y + d x V u x = -£Tq^(£ + 2 7) 



du x du y 

dx dy 

du x duy 

dx dy 



where 



Odd 
at ox ay 



d 2 

dx 2 dy 2 



(2.5a) 
(2.5b) 

(2.6) 



Finally, note that u z is give n by the components u x and u v through the continuity equation. 

The temporal mode ( [2.3] ) may be decomposed upon a basis of such modified Fourier modes (2.4). Let us expand 
the spatial part (v w , q u ) of (2.2) in the usual Fourier modes along the z direction : 



(vo,,gu. 



(v(a;,y,/c ), 9 (x,y,A:o))e l ' Co ^fco . 



(2.7) 



This expansion can be viewed as a superposition at time t — of generalized Fourier modes provided the following 
initial condition is satisfied : 



(u(x,y,0,k o ),p(x,y,0,k o )) = (v(x,y,k ),p(x,y,k Q )) 



(2.8) 



Accord ing to (2.4), each generalized Fourier mode evolves independently: an alternative expression for the temporal 
mode (2.3) is thus provided for all t 



(u w ,p a 



(u(x, y, t, k ),p(x, y, t, feo )) e * e -^, e -^(i-e— )/ 2o _ 



(2.9) 



The above expression is consistent with (2.3) if the following equality holds at any x and y locations, time t and 
wavenumber kg : 

(u(x, y, t, k ),p(x, y, t, ^K^ 1 - 6 "^ 27 = (v(as, y, fc e" 7 *), q(x, y, k^)^ e~ wt . (2.10) 

In the following section, equality ( p. 10 ) is shown to be valid only for fco = 0: discrete temporal modes are bound to 
be two-dimensional. 
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III. THE Z DEPENDENCE OF A THREE-DIMENSIONAL TEMPORAL MODE 



First a cut-off wavermmber k c above which v(x, y, ko) and g(x, y, ko) vanish, should exist. In deed, assume that such 
a cut-off does not appear, large wavenumbers are the n pre sent in the spatial spectrum of (2.3). It is thus possible to 
take the simultaneous limits hot large and t small in ( 2.10 ). The right-hand side of this equation becomes 



(v(x,y,fc e 7 *), q(x, y, k e 7 *))e 



(v(x,y,k ),q(x,y,k )) 



(3.1) 



In order to estimate the left-hand side of ( 2.10| ), the behavior of u(x, y, t, ko),p(x, y, t, ko)) is evaluated using equations 
(2.5a,b). Two cases are to be considered according to the characteristic spatial variations of u x and u y in the x and y 
directions. When, for la rge axial wavenumber ko, these components evolve over spatial scales independent of fco, the 
right-hand side of (2.5a,b) can be neglected and the lead ing order time evolution is independent on fc . This means 
that, for kgt large and t small, the left-hand side of ( 2.10| ) reads 



(u(x,y,t,k ),p(x,y,t,k ))e vk ° {1 e 27 ' )/27 ~ (u^x, y, 0),Poo(x, y, 0))e 



(3.2) 



On the contrar y, wh en u x , u y evolve over spatial scales comparable to 1/fco, the correct expression is found when 
introducing in ( [2.54 b) new fast variables 



k x ; y = k y ; t = k Q t, 



(3.3) 



in addition to x and y, and thereafter expanding the solution with respect to 1/fco- At leading order, homogeneous 
equations are obtained : 



dux I 

-w-"\ 


' d 2 


d 2 N 






du x 


f d 2 


d 2 


[dx 2 ' 


h dy 2 
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dx 2 
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du x du h 
dx 
du x 
dx 



du y 
dy 



(3.4a) 
(3.4b) 



For t small, the term e 27 ' may be taken equal to 1 and the general solution of ( 3.4a ) be written in the fast variables 
x and y via the usual Fourier decomposition 



(u x ,u y ) ■■ 
For fixed x and y and fco^ 



(4°) (k x ,k y , x, y), (k x ,k y , x, y))e-^ k2 * +k2 y*e tk * x+tk yydk x dk y . 
•oo, 



the above expression (3.5) is evaluated by the steepest descent method 

_ v 47T . fnl ,„ „ \ ^m\,~ „ ^^ 



vk 2 t 



(u x o \0,Q,x,y),u y °\0,O,x,y)) 



Similar behaviors can be obtained for u z and p. For this case, the right-hand side of ( |2.10 ) now reads 



(u(x, y, t, k ),p(x, y, t, k^e-^ 1 '^/^ ~ (0, 0, x, y),pW ( , 0, x, y)) 



47T 

vk^t 



-uk n t 



(3.5) 



(3.6) 



(3.7) 



Introducing into equation (2.10), estimates ( p.l[ ) and ( |3.2| ) or (3.7) leads to an inconsistency for large wavenumbers 
ko. This contradiction implies that the spatial spectrum of ( [2.3] ) vanishes for sufficiently large wavenumbers : there 
hence exists a cut-off wavenumber k c such that (v(x, y, ko), q(x, y, ko)) = for fco > k c . 

Consider now a wavenumber < fci < k c such that (v(x, y, ki), q(x, y, k\)) ^ (0,0) and a time t\ such that 
ko = /cie 7 * 1 > k c . Equation ( 2.10| ) then implies that (ii(x, y,t%, ko),p(x, y, t\, ko)) which leads again to a 
contradiction. Indeed, the initial condition (u(x, y, 0, ko),p(x, y, 0, fco)) is equal to zero since fcp > k c , the quantity 
(u(x,y,t,ko),p(x,y,t,ko)) thus remains zero for all t > since it is governed by equation ( [2.5a[ b). This imposes 
that (v(x, y,ki), q(x, y, fci)) = (0, 0) for all fci / 0. Discrete temporal modes are then independent upon z i.e. purely 
two-dimensional ones. For the axisymmetric Burgers vortex, those have been computed by Robinson & Saffman 
and Prochazka and Pullin These authors showed that these modes are damped. 

Needless to say, the general stability problem for Burgers vortex or any stretched vortex, cannot be solved only 
looking at the discrete part of the temporal spectrum. In particular, the ana lysi s of the continuous spectrum need to 
be considered. However, the process of "bidimensionalization" displayed in ( |2.4| ) can be used again in that case. 
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